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Abstract. A spread code is a set of vector spaces of a fixed dimension over a 
finite field ¥ q with certain properties used for random network coding. It can 
be constructed in different ways which lead to different decoding algorithms. 
In this work we present a new representation of spread codes with a minimum 
distance decoding algorithm which is efficient when the codewords, the received 
space and the error space have small dimension. 



1. Introduction 

In network coding one is looking at the transmission of information through 
a directed graph with possibly several senders and several receivers pQ. One can 
increase the throughput by linearly combining the information vectors at interme- 
diate nodes of the network. If the underlying topology of the network is unknown 
we speak about random linear network coding. Since linear spaces are invariant 
under linear combinations, they are what is needed as codewords [7J. It is helpful 
(e.g. for decoding) to constrain oneself to subspaces of a fixed dimension, in which 
case we talk about constant dimension codes. 

One class of constant dimension codes is the one of spread codes. These codes 
have maximal minimum distance and are optimal in the sense that they achieve the 
Singleton-like bound on the cardinality of network codes. They can be constructed 
with the help of companion matrices of irreducible polynomials, as explained in 

In this work we translate this construction to an extension field setting and 
evolve a decoding algorithm for spread codes in this setting. The complexity of this 
new algorithm depends on different parameters than the algorithms of O [7j lllj , 
which are also applicable to spread codes. Therefore, depending on the network 
setting applied on, the new algorithm has an improved performance. 

The paper is structured as follows: In Section [5] we give some preliminaries 
on random network coding and constant dimension codes. The main results of 
this work are found in Section^ where we first show how to translate the spread 
code construction into a different setting and then explain how decoding can be 
done in this setting. We study the complexity of the decoding algorithm and give 
comparison to other known decoding algorithms in Section 01 
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2. Preliminaries 



Let F g be the finite field with q elements, where q is a prime power. We denote 
the set of all subspaces of F™ by PG(n, q) and the set of all fc-dimensional subspaces 
of F™, called the Grassmannian, by Q q (k,n). The general linear group GL n is the 
set of all invertible n x n-matrices with entries in ¥ q . Moreover, the set of all 
k x n-matrices over F g is denoted by Matkxn- 

Let U G Matkxn be a matrix of rank k and 

U = rs(Z7) := row space(C/) £ Q q {k,n). 

One can notice that the row space is invariant under GL^-multiplication from the 
left, i.e. for any T £ GL k 

U = vs(U) = vs(TU). 

Thus, there are several matrices that represent a given subspace. A unique repre- 
sentative of these matrices is the one in reduced row echelon form. Any k x n-matrix 
can be transformed into reduced row echelon form byaTs GLk- 

A subspace code is simply a subset of PG(n,q) and a constant dimension code 
is a subset of the Grassmannian Q q {k,n). 

The subspace distance, given by 

d s (U, V) = dim(W + V) - dim(W n V) 

for U, V two subspaces of F™, is a metric function on PG(n, q). It induces a metric 
on Q q (k, n) by 

d s (U, V) =2k - 2 dim(W n V) 

for any U,V £ Q q (k,n). The minimum distance of a subspace code C C PG(n,q) 
is defined as 

d{C) := mm,{d s (pi, V) \ U,V G C, U V}. 

The subspace distance is a suitable distance for coding over the operator channel 
[7] , where errors and erasures can be corrected. An error corresponds to an inserted 
erroneous vector, i.e. an increase in dimension, whereas an erasure is a decrease 
in dimension of the code word. The error-and-erasure capability of a code C C 
PG(n, q) with minimum distance d(C) is 

_ d(C) - 1 
L 2 

Different constructions for constant dimension codes can be found e.g. in [31 [6| \7\ 

Eng. 

In the case that k divides n one can construct codes with minimum distance 2k 
and cardinality q qk Z\ ; called spread codes [10] . One construction for spread codes 

is the following. Let p(x) — ^2^ =0 PiX l £ F g [x] be an irreducible monic polynomial 
of degree k. Its companion matrix P £ Math 
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Then the F g -algebra of such a companion matrix ¥ q [P] is a finite field and the set 



is a spread code. Since we want to work with the reduced row echelon form as a 
unique matrix representation of the vector spaces we assume that the first non-zero 
block from the left is the identity matrix. These codes are optimal because they 
achieve the Singleton-like bound [7] on the cardinality of constant dimension codes 
for a given minimum distance. 



3.1. Translation of the Construction. We now translate the construction 
of spread codes of [101 from the companion matrix to an extension field setting. 
Since they exist for any degree over ¥ q , we choose primitive polynomials and their 
companion matrices for the spread code constructions. 

For the remain of this paper assume that k\n and let I = n/k. Moreover, let 
a G ¥ q k be a primitive element of ¥ q k and /3 G ¥ q n a primitive element of ¥ q n 
as an extension field of ¥ q k. The polynomial p(x) G denotes the primitive 

polynomial of degree k such that p{a) = and P G its companion matrix. It 
holds that ord(P) = 9 fc -l| Lemma 6.26]. 

Denote by : F* — > ¥ q k and (f>^ : ¥ l qk — > ¥ q n the standard vector space 
isomorphisms: 



{rs [ B B x 



B f -x ] eg q (k,n) | B t G ¥ q [P] } 



3. Spread Codes in Extension Field Representation 




i=0 



Vl-l) 



for (u , ■ ■ ■ ,itfc-i) G ¥ k q and (v , . . . G F l qk . 



Proposition 1. Let e\, . . . , e n be the standard basis o/F™. Then 



i-i 



is a basis of ¥ q n over ¥ q and 




is a vector space isomorphism. 



Proof. Define :W™ —>¥ l q 



,k i 



4> {k) (ui, . . . , v n ) := U (fe) (vi,..., v k ), . . . , 4> {k) K-fc+i, . . . , v n )\ 
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Then </>W , (j)( k ) are vector space isomorphisms and cf> = <j>^ o </>( fe ) satisfies the fol- 
lowing diagram 

F£ * >F 9 n . 





q K 

□ 

Note, that <j> can be applied on sets of vectors (e.g. vector spaces) element-wise. 

Lemma 2. Denote by P[i] the i-th row vector of P. Then 

<j>W{P h \i]) = a h+i - 1 

for i = 1, . . . , k and h = 1, . . . , q k — 1. 

PROOF. It is easy to see that (j)^ k \P[i\) = a % for i G {1, . . . , fc}. Moreover, (j)^ 
is commutative with the multiplication with P and a: 

k 

4>^{uP)=Y J UiPij^- 1 
»=i y =1 

k 
i=l 

for all u = (wi, . . . ,Ufe) G F*. 

=► (fc) (P h [*]) = ( /.( fe )(PHP /l - 1 ) = 0( fc )(P[*])a fc - 1 = a h+l - x 

□ 

Recall that the spread code elements are of the type 

rs[ B B 1 ... Bi-! ] eg q {k,n) 

where each block Bi is an element of ¥ q [P] and the first non-zero block from the 
left is the identity. 

Theorem 3. Define 

f if B 1 = 

ir-= [ a h i{Bj = P h 

Then 

i-i 

0(rs[ Bo B 1 ... Bi-! ])=F,» .^7^'. 

3=0 

Hence, we can uniquely identify each spread code element by the respective 7 = 
(70,..., 71-1) e^- 
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Proof. Denote by Bj[i] the i-th row vector of the block Bj. From Lemma [2] 
we know that (jP^ (Bj [i]) = a 4_1 7j. The power of j3 corresponds to the position of 
the block Bj , thus in general <p maps the i-th row of the whole matrix to 

2-1 Z-l 

]T a*- 1 ^ = a*- 1 IjP 3 Vi = 1, . . . , k. 

j=0 j=0 

As a is a primitive element of ¥ q k , the elements of the vector space are exactly 
mapped to ¥ q k ■ Y^jJoljft- ^ 

Corollary 4. A spread code C C Q q (k, n) constructed as before is isomorphic to 

MM)- 

Proof. Since a spread code covers the whole space and <j> maps a code word 
to an F„*-linear subspace with basis vector 7, the statement holds. □ 

Spreads of this type are also known as ¥ q -linear representations of PG(l,q k ) 
[2] or Desarguesian (k — \)-spreads [8] p. 12]. 

3.2. Decoding. One can now use this structure for the decoding procedure 
in this representation. 

First assume only erasures and no errors happened during transmission. Then 
any received vector space TZ with dim(7?.) > 1 can be decoded to its closest code 
word, since the number of errors and erasures is less than or equal to k — 1 = 
2fc ~ 2 = t. For decoding choose an element of the received space r £ TZ and compute 
7 = (70, . . . , 7/_i) S ¥ l qk such that 

I-i 
3=0 

for some i. For this, divide r into I blocks of size fc, r%,. . . ,r;, and find the first 
non-zero block, denoted by r s . It holds that r s — (j)^ (a 1-1 ), since the first non- 
zero block is the identity matrix in the construction. Then 7 can be computed by 
I divisions in ¥ q k of the image under </>W of each block by r s . 

7 = (0 w (n)0«(r.) _1 . ■ ■ ■,0 (fc) (n)* (fc) (r.)- 1 ) 

Note, that one does not need to compute the discrete logarithm to find the i. It is 
enough to compute 7 to identify the code word. 

Example 5. In the [3, 6]-spread code over F2 let r = (1 10 1 101) be a received vector. 
It holds that <p(r) — 1 + a + (3 + a 2 (3. The first three entries tell you that you have 
to divide by 1 + a to compute 7: 

(1 + a + $ + a 2 f3)(l + ay 1 = 1 + (1 + a)/3 

Hence, 7 = (1, 1 + a) = (1, a 3 ), which identifies the codeword 

rs [ I I P 3 ] = rs 
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But what if errors were inserted? Let U £ C be the sent code word and denote 
by k' the dimension of the received vector space 1Z £ PG(n, q). For correct decoding 
it has to hold that 

<S=> k + k' -2 dim(W n K) < k - 1 
<=*> dim{Ur\TZ)>^-, 

therefore one needs to find f^-g^] linearly independent elements of 1Z with the 
same respective 7, called 7 max - Then we decode to the codeword 

l— 1 

3=0 

Since we do not know if any or which of the elements of 1Z are erroneous, one needs 
to examine not only a basis but all elements of 1Z. 

A first basic decoding algorithm in this extension field representation is given 
in Algorithm [T] All field operations are done over ¥ q k . 

Algorithm 1 Basic decoding algorithm. 

Require: the received vector space 7Z £ PG(n,q), k' — dim(7£) 
for each v £ 1Z do 

divide v into blocks vq, . . . , vi~i of length k 

v s := the first block from the left with non-zero entries 

a:=(0W(t,,))-i 

store 7„ := (ft® (v ) -a,..., tf>W • a) 

end for 

7max := the element of highest multiplicity in {^ v \v £ 7Z} 

if there are > linearly independent v £ 1Z such that j v = 7 max then 

return -1 (F,* • X^-=o 7maxj/3 J ') 
else 

return "not decodable" 
end if 



3.3. Improvements on the Algorithm. We improve the algorithm by sys- 
tematically choosing the linear combinations of the basis vectors of the received 
space to work with. For it, note that errors are canceled out in some linear combi- 
nations of elements, as illustrated in the following example. 

Example 6. Assume U £ C was sent and consider two elements of the received 
space r\,T2 £ 7Z containing the same error e £ F™, i.e. 

ri=/ J A u u + e , r 2 = ^ Mu" + e. 

ueu ueu 

Then 

n + (q- l)r 2 = ^2 X u u + e-^2 \« - e = (A„ - n u )u £ U. 
ueu ueu ueu 

Let us generalize this idea to arbitrary numbers of errors. 
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Proposition 7. Let ui, . . . , itfc 6 F™ &e a &aszs o/ i/ie seni code word U £ <7 9 (fc, n) 
and ri,...,rv G F^ a ftasis o/ f/ie received space 1Z. Assume / < fc' linearly 
independent error vectors were inserted during transmission, i.e. TZ — W ' ®£, where 
W is a vector subspace of IA and £ is the vector space of dimension f spanned by 
the error vectors. Then the set 

|5^Vi|A i €F„/c{l,...,fc , } ) |J| = /}| 

contains fc' — f linearly independent elements of Li. 

Proof. Inductively on /: 

(1) If / = 0, then n,..., tv eU. 

(2) If / = 1, assume r\, . . . , rj ^ W and n+i, ■ • • , ?V S W ■ Then there exist 
\ e¥ q , m e¥ q \ {0} such that 

fc 

r i = '^2 hjUj + Mi e V i = 1, . . . , / 

2=1 

where e € £ denotes the error vector. Hence \f i,h = I, . . . ,1 

k 

n + r h = y^(Ajj + Ahj)uj + (/Xj + Hh)e 

2=1 

fe 

=> r, + (-^i/i^ 1 )^ = y^(A»3 + A/, 3 -)itj e W. 

2=1 

Then the elements n+i, . . . ,?v,ri + (-^i^ 1 )r 2 , . . . ,n + (— /ii/^~ )r/ are 
k' — 1 linearly independent elements without errors. 

(3) If more errors, say ei,...,e/, were inserted, then one can inductively 
"erase" / — 1 errors in the linear combinations of at most / elements. 
Write the received elements as 

k f 

n = x a u 3 + Mi2 e 2 V i = l, . . . , fc' 

2=1 2=1 

with Xij,fJ,ij G F g . Assume /«i/, . . . ,fXif ^ and ■ • ■ i /"fe'/ = 0, i.e. 

the first Z elements involve e/ and the others do not. 

From above we know that the linear combinations of any two el- 
ements of r\ , . . . , r\ include I — 1 linearly independent elements with- 
out ef. Denote them by mi, . . . , roj—i. Naturally these elements are 
also linearly independent from r; + i, . . . , rv. Use the induction step on 
mi, ... , m;_i, n+i, . . . ,rv to get fe' - 1 - (/ - 1) = fc' - / linearly inde- 
pendent elements without errors. 

□ 



Corollary 8. In the setting of Proposition^ assume ds{1Z,U) < t, i.e. 1Z is 
decodable. Then there are at least \ k j" 1 ] linearly independent elements ofbt in the 

C:=\j2^n I *i GF ? ,/C {!,..., k'},\I\ = f 



2 

set 



iei 



s 
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Proof. Let / denote the number of erasures. Then / = / + k — k' and thus 

k' - 1 



f + f <t 2f + k-k'<k-l / < 

With Proposition [7] it follows that C contains k' — f > linearly independent 
vectors of the sent vector space. □ 



We use this fact to modify Algorithm [T] as follows: We choose a basis 7*1, . . . , ry 
of the received space 1Z € PG(n,q) and compute r y ri for i = l,...,fc'. Then we 
compute the respective 7 of all linear combinations of two basis elements, then of 
three elements etc. As before we can stop the process and decode to a code word as 
soon as we have more than or equal to \ k j" 1 ] linearly independent elements with 
the same 7. This way, if / errors occurred, we do not have to consider all elements 
of 1Z but only the linear combinations of at most / of the basis vectors. 

Moreover note, that a linear combination of elements with the same 7 is always 
another element with 7. Since we need to find linearly independent elements, it is 
therefore enough to check only combinations of elements with different respective 
7's. 



4. Complexity 

For a better understanding of the complexity of the algorithm we first consider 
binary spread codes and then generalize it. Note, that the algorithm works for 
received spaces of arbitrary dimension. 

If k' is the dimension of the received space and / is the dimension of the 
error space £, the algorithm computes the sums of at most / + 1 basis vectors, 
which are (j k + i) many. For each sum it proceeds with an inversion and at most 
x — 1 multiplications over F 2 fc . The complexity of inverting is upper-bounded by 
0{k 2 ) over F 2 and the one of multiplying by O(fclogfc) over F 2 using the FFT [H 
Chapter 8.2]. Using the approximation (j k +1 ) ~ (j+iji 1 the overall complexity is 
upper- bounded by 0{nkk'f +1 ) over IF 2 ■ 

Over F q one needs to consider not only sums but F g -linear combinations. Thus 
we get an upper bound of (j+jj combinations to check. Hence, the overall complex- 
ity is upper-bounded by 0(nk(qk')f +1 ) over F 9 . 

The complexity reduces when some of the generators of the sent codeword are 
not influenced by the errors since in this case the algorithm has to check only linear 
combinations of a smaller amount of basis vectors of the received space. 

In the following we compare this complexity with the one of the spread decoding 
algorithm shown in [5] and the decoding algorithms for Reed-Solomon like codes 
contained in [7j and in the case of q = 2 and k = k' . In [5] the authors 
present a minimum distance decoder for the their spread code construction. The 
complexity of their algorithm is 0((n — k)k 3 ). If the dimension of the error space 
is minimal the two algorithms perform the same. When applied to spread codes 
the complexities of the algorithms presented in [7] and are 0(n 2 (n — k) 2 ) and 
0(k(n — fc) 3 ), respectively. The algorithm proposed in this work performs better 
if the dimension of the codewords, of the received space and of the error space are 
small. 



SPREAD DECODING IN EXTENSION FIELDS 



9 



5. Conclusions 

In this work we consider a new construction for spread codes, which are codes 
with optimal cardinality and error correction capability. The construction is based 
on the representation of code words via a unique element 7 £ ¥™ k . We present a 
minimum distance decoding algorithm which works by finding this unique 7. The 
performance of the algorithm is mainly based on two properties: the operations 
are done over ¥ q k instead of ¥ q ™ and the elements tested for finding the 7 are only 
the linear combinations of / of the basis vectors of the received space instead of all 
elements of the vector space. As a result we obtain a decoding algorithm with a 
good performance when the dimension of the codewords, of the received space and 
of the error space are small. 
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